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Abstract 

We consider inviscid limits to shocks for viscous scalar conservation 
laws in one space dimension, with strict convex fluxes. We show that we 
can obtain sharp estimates in L 2 , for a class of large perturbations. Those 
perturbations can be chosen big enough to destroy the viscous layer. This 
shows that the fast convergence to the shock does not depend on the fine 
structure of the viscous layers. This is the first application of the relative 
entropy method developed in [22] , [23] to the study of an asymptotic limit 
to a shock. 
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1 Introduction and the main result 

For any strictly convex flux function A E C 2 (R), we consider the family of 
viscous scalar conservation laws in one space dimension: 

d t U + d x A{U) = ed 2 xx U for (>fl,iel, 
£7(0,2;) = U {x) for i£l, 

for any e > and Uq € L°°. Global unique solutions to (1) have been con- 
structed by Hopf [17] and Oleinik [28]. The inviscid case, e = 0, is covered 
by the theory of Kruzkov [20]. Kuznetsov showed in [21] that, for fixed initial 
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values Uo, the solutions of (1) converges, when e goes to zero, to the solution of 
the inviscid scalar conservation law (the equation (1) with e = 0). He showed 
also that, for general initial data, the optimal rate of convergence in L 1 is y/e. 
In the case of the convergence to a shock, however, the rate is better. This is, 
usually, linked to the formation of layers. 

In this paper we consider the asymptotic limit for general initial values. We 
are particularly interested in the cases where the initial values carry too much 
entropy for the structure of the layer to be preserved asymptotically. The shocks 
solutions of the inviscid case (e = 0) can be described as follows. Consider two 
constants Cl > Cr, and the associated function defined by 

r, / x J C L if x < 0, 

So(x) = < (2) 

! Cr it x > o. 

Then, the Rankinc-Hugoniot conditions ensures that the function 

Qf ,x A{Cl) - A(C R ) 

S {x - at) with a := , (3) 

is a solution to the inviscid equation (1) with e = 0. The condition Cl > Cr 
implies that they verify the entropy conditions, that is: 

dtn{U) + d x G{U) < 0, t > 0, x e R, 

for any convex functions 77, and G' = rj'A'. 

Our main result is the following. 
Theorem 1.1. Let C L > C R and U a € L°°(R) n BVi oc (R) be such that 

(U -S a )eL 2 (R) and (^C/ )+ £ L 2 (I). 

Then, for any T > 0, there exist £0 > and C* > such that the followings 
hold true. 

I. For any U solution to (1) with < e < Sq, there exists a curve X £ 
i°°(0, T) such that X(0) = and for any < t < T: 

\W(t) - S(t)||£ a(R) <\\U - S \\h m + C* £ log(l/ £ ), (4) 

where S(t,x) :— So(x — X{t)), and So is defined by (2). 
II. Moreover, this curve satisfies 

\X(t)\<C* and (5) 

\X(t) at\ 2 < C*t 2 '* (||C/ - Sol| 2 L2(R) + elog(l/e)) . (6) 
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777. The constant e$ depends only on ||(^£^o)+||l 2 ; Cl, Cr ; \\Uq\\l°° and the 
flux function A, while C* depends only on the same set as well as T. 

Remark 1.1. Note that our estimates do not depend on any local BV norms of 
Uq. The assumption Uq e BVi oc ensures that -^Uq is a Radon measure. Hence, 
(^f/o)+ is also a Radon measure, and the condition (^{/o)+ € L 2 makes sense. 

This result shows a rate of convergence slightly worse than e (to the log), for 
the inviscid limit to a shock, measured via the L? norm (squared). In the case of 
the limit to a regular solution of the inviscid case, the rate of convergence is y/e 
(see [32], for instance). We also refer to [21], Goodman and Xin [16], Bressan, 
Liu and Yang [6] , Lewicka [24] , Bressan and Yang [5] , Christoforou and Trivisa 
[10]. 

An easy layer study shows that e is the optimal rate for shocks with special 
initial data. Indeed, one can construct an associated steady viscous layer (see 
for example II' in and Oleinik [18]) Si solution to 

J A(St) - A(C L ) - a{St - C L ) = S[, x e R, 

I lim Si = C L , lim Si = C R . ^' 

k x^? — oo x— 5-+00 

It is easy to show that S\{{x — crt)/e) is a solution to (1) with initial data 
S\{x/e). In this case, the rate of convergence is of order e since: 

/ \S 1 ((x-at)/e)-S (x-at)\ 2 dx<e [ |Si(x) - S (a;)| 2 dx = Ce. 

This layer study can be extended to the case of small initial perturbation where: 

f \U (x) - S {x)\ p dx < Ce, 

for a 1 < p < oo. In this case, for a solution U to (1), we can consider 

V(t,x) = U(et,sx), 

and study the asymptotic for large time. The function V is a solution to the 
equation 

d t v + d x A(v)-d 2 xx v = o, 

V(0,x) = U(0,x/e). 

The convergence to Si, up to a (constant) drift, in this setting, has been ex- 
tensively studied (see for instance [18], Freistiihler and Serrc [14], Kenig and 
Merle [19]). In this situation of small perturbation of the initial shock, those 
results show that the convergence with rate e for the system (1) is due to the 
asymptotic limit in large time of the layer function U{-/e) to Si((- — at)/e). 
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This layer study, however, collapses when 



f \U (x) - S (x)\ 2 dx 



> e. 



In this situation, there is too much entropy for the asymptotic limit of the layer 
structure to be true. The physical layer may be destroyed. Theorem 1.1 shows 
that, nevertheless, the sharp convergence (up to the log) still holds. 

Taking a limit as e goes to in Theorem 1.1, we recover the L 2 stability 
of shocks (up to a drift) first showed by Leger in [22]. Note that the stability 
result has to be up to a drift which depends on the solution itself (and may be 
not unique). This feature is also true for our result. The drift cannot be taken 
constant, as in the case of the layer problem. 

Our result is based on the relative entropy method first used by Dafermos 
and DiPerna to show L 2 stability and uniqueness of Lipschitzian solutions to 
conservation laws [11, 12, 13]. They showed, in particular, that if U is a Lip- 
schitzian solution of a suitable conservation law on a lapse of time [0,T], then 
for any bounded weak entropic solution U it holds: 



for a constant C depending on U and T. 

The relative entropy method is also an important tool in the study of asymp- 
totic limits. The main idea is that convergence holds thanks to the strong sta- 
bility of the solutions of the limit equations. Roughly speaking, if we have good 
consistency of e models, with respect to the limit one, then non linearities are 
driven by the strong stability of the solution of the limit equation. Applica- 
tions of the relative entropy method in this context began with the work of 
Yau [33] and have been studied by many others. For incompressible limits, see 
Bardos, Golse, Levermore [1, 2], Lions and Masmoudi [25], Saint Raymond et 
al. [15, 30, 26, 29]. For compressible models, see Tzavaras [31] in the context 
of relaxation and [4, 3, 27] in the context of hydrodynamical limits. However, 
in all those cases, the method works as long as the limit solution is Lipschitz. 
This is due to the fact that strong stability as (8) is not true when U has a 
discontinuity. It has been proven in [22, 23], however, that some shocks are 
strongly stable up to a shift (see also related works from Chen and Frid [7, 8] 
and Chen, Frid and Li [9]). This article is the first extension of those results 
of stability, to the study of asymptotic limits to a shock. This is a part of the 
program initiated in [32] . 

The result can be extended to any entropy in the following way. Fix any 
strictly convex function r\ e C 2 as an entropy. We define the associated relative 
entropy functional as 





rj(x\y) := r)(x) - r](y) - r)'(y)(x - y). 
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We then have the following extension. 

Theorem 1.2. Consider a strictly convex entropy functional rj e C 2 (M). Let 
C L > C R and U G L°°(R) n W ioc (R) oe sucft tfiat 

(C/ -5 )e£ 2 (R) and (^o)+ e L 2 (K). 

Then, for every T > 0, there exist Eq > and C* > smc/i f/iaf i/ie followings 
hold true. 

I. For any U solution to (1) with < e < So, there exists a curve X G 
L°°(0,T) smc/i t/iat X(0) = 0, and for any < t < T, and for any a 
verifying e < a < Eq, we have: 



r]{U{t,x)\S{t,x))dx< [ 7]{Uo(x)\S {x)) dx + C* 

{\x-X(t)\>C*a} JR 

(9) 



where S(t,x) := Sq(x — X(t)), and Sq is defined by (2). 
II. Moreover, this curve satisfies 

|A"(t)|<C* and (10) 



\X(t)-at\ 2 < C*t 2/3 ( f n{U (x)\S (x))dx 
\Jr 



+ elog(l/e) . (11) 



The constant So depends only on ||(^-[/o)+||l 2 ; Cl, Cr, ||^o||loo 7 the flux 
function A and the entropy functional r\, while C* depends only on the 
same set as well as T. 

Theorem 1.1 is a direct application of Theorem 1.2 with n(x) := x 2 , and 
a = elog(l/e). Indeed, in this case we have r](x\y) = (x — y) 2 , and 

/ r]{U(t, x)\S(t, x)) dx < C\{\x - X(t)\ < C*a}\ < CC*a. 

J{\x~X(t)\<C"a} 



For the rest of the paper, we will assume that the initial value U lies not only 
BVi oc but also C x . It allows us to work with smooth solutions U € C 1 ([0, T] xR). 
The general BVi oc case can be obtained by a density argument. 

The idea of the proof is to study the evolution of the relative entropy of the 
solution with respect to the shock, outside of a small region centered at X(t) 
(this small region corresponds to the layer localization): 



pX(t) — Se p 00 

I r){U{t,x)\C L )dx+ / n{U{t,x)\C R )dx. (12) 

-00 JX(t)+8e 
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The change in time involves two effects. One is due to the hyperbolic part of 
the equation, and the second involves the parabolic part (or order e). In [22], 
it was shown that, for the hyperbolic case e = 0, with zero layer width 5 = 0, 
this quantity is non-increasing when we choose wisely the drift X(t). When 
considering the viscous term, the layer with width Se is introduced to mollify 
its effect. The idea is then that the stability induced by the hyperbolic part is 
enough to compensate the parabolic one, provided that we consider a layer fat 
enough. For technical considerations, we will hx S = log(l/e). The drift X(t) 
is still chosen with respect to the hyperbolic part of the equation in a similar 
way as in [22]. The stability effect is preserved, despite the non zero layer 
width, thanks to a monotonicity property induced in the layer by the additional 
assumption (^f7 )+ € L 2 (R). 

2 Evolution of the relative entropy 

For S > 0, we consider a Lipschitz nondecreasing function (j>$ to localize the 
layer, verifying 



4>s(x) 




it x < 0, 
if x > S. 



To get the optimal result, we will later fix a special function (see (19)). 

For any fixed 5 > and X £ C 1 ([0,T]), we are interesting in the evolution 

of 

/OO 
${\x - X(t)\/e) V (U(t, x)\S(t, x)) dx, (13) 
-oo 

where the function S depends on X as in the statement of the theorems. A 
special value of 5 (depending on e), and of the function X will be chosen later. 
Note that H(t) is controlling the quantity (12). From now on, we drop the 
subindex 5 for <j>s for a simple presentation. 

Let us denote F(-, ■) the flux of the relative entropy rj(-\-) defined by 

F(x, y) := G(x) - G(y) - r/'(y)(A(x) - A{y)). (14) 

Note that 

d uV (U\C)=r 1 '(U)-r ] '(C), 

d v F(U, C) = G'{U) - v'(C)A'(U) = ( V '(U) - r/(C))A'(U). 

So, for any solution U of (1) and any constant C, multiplying (1) by rj'(U) — 
rf{C), we get 

d tV (U\C) + d x F(U, C) - e{rf{U) - v'(C))d 2 xx U (15) 
We have the following lemma. 
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Lemma 2.1. T/ie function H, defined in (13), satisfies the following on (0,T) 



dt 



(*) = 



JX(t)-Se £ v £ 7 v £ 7 L 



X(i)-«e fc N e 7 V 



da; 



_ f X(t)+S£ H (£_J[Wy [x(t)r,(C/(t,x)|C«) - F(I/(t,s),C fl ) 



X(t) 



.x - X(t) 



d 2 xx U(t,x)( V '(U(t,x)) - rj'(C R ))dx 



'■- ( L )Hyp + ( L )Dif + ( R )Hyp + ( R ) Dif 
Proof. First we split the above integral into the two parts: 



H(t) 



X(t)\ 



n{U\S)dx 



(W =£d i m )]^W £z r m )]>^ 5 ) <fa 



x-x{ty 



'x-x(ty 



:= H L + H H . 

To compute ^{H L ), we put C = Cl in (15), multiply by 
and integrate in x. Then we have 



-x+X(t) 



-x + X(t) 



£< fll >-/>(K . 



+ /_ « 

+ £ 



oc 

r L( -* + A w ) 1 *<%M4(u) - v '(c L ))dx 

J — oo £ 

f x "> (?) . ,(z£±£« y(zE±£») [t mum - F{ u, c L) 

Jx(t)-Se £ V £ 7 V e 7 L 



'-X + X(t)\l2, 



J7(I7|Ci)da; 
F(U,C L )dx 



dx 



x(t) \,(~ x + x ( t Y 



r w r / 

7 K 



i)] d 2 xx U(n'(U) - n\C L ))dx = (L) Hyp + (L) mf . 
We get the result for f t (H R ) = (R) Hyp + (R) Di f in the same way. □ 
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3 Study of the contribution of the hyperbolic 
structure 

In this section, we show that by choosing a special drift function X(-), the hy- 
perbolic effects become nonpositive. This will be used in section 4 to control 
the parabolic effects. 

Following [22], we define the normalized relative entropy flux /(•,•) by 

,/ x F(x,y) 
r){x\y) 

We have the following properties. 

Lemma 3.1. For any L > 0, there exists a constant A > 7 such that for any 
x,y with \x\, \y \ < L, we have 



1/A < r]"(x) < A, 

-L(x-y) 2 < V (x\y) <^A(x-y) 2 , 

\F(x,y)\<A(x-y) 2 , 
< (d x f){x,y) < A, 
l/A<(d y f)(x,y). 

The proof of this lemma can be found in [22] . 

We now define the shift function X. It is the solution of the following O.D.E. 
\ X(0)=0 

Note that for any e > 0, U £ C^^T] x M) (since U € C^R)). The existence 
and uniqueness of X comes from the Cauchy-Lipschitz theorem. 

First, X is Lipschitz, since we have from Lemma 3.1 



F(u(t,X(t)) 


C L +C R \ 
2 J 


v(u(t,x(t)) 


C L +C R \ 

2 ) 



X{t)\ < — ) f- < 2A 2 (18) 



where we used the fact ||{7(£)||l°o < ||C/o||l°° for t > 0. It proves (10). 

Note that from the definition of X, if the the velocity U were constant in 
the layer (that is U(t,x) <~ U(t,X(t)) for — Se < x — X(t) < <5e), then, from the 
last property of Lemma 3.1, we would have directly that 

(L) Hyp + (R) Hyp < - ° L a CR (r,(U(t,X(t))\C L ) + V (U(t,X(t))\C R )). 
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However, this is too much to hope, since the layer characterize the region where 
the function U(t, •) is expected to drop from Cl to Cr. We still can show that 
the hyperbolic terms are negative, provided that the behavior of U in the layer 
is not too much oscillatory (the values can drop, but not much bounce back). 
This last property of U is proved in the following lemma which can be seen as 
a weak version of the Oleinik's principle . 

Lemma 3.2. \\(d x U(t)) + \\ L 2 {M) < \\(-£Uo)+\\^(R) for any t > 0. 

Proof. We differentiate (1) w.r.t. x, multiply (d x U)+ and integrate in x to get 







d t d x U + A"(U)\d x U\ 2 + A'(U)d 2 xx U - sdl xx U 



dx 



-d t ([(d x U) + ] 2 )+A"(U)(d x Uf + 

+ A'(U)d x ( [{dxU 2 )+]2 ) + e\d x ((d x U) + f 



dx. 



Then, we use the integration by parts to get 



/ 



1 



-d t ([(d x U) + } 2 ) + A"(U)(d x U)l+e\d x ((d x U). 



dx 



> 



1 d 
Ydt 



J [{d x U)+] 2 dx. 



□ 



We now prove the main proposition of this section. 



Proposition 3.3. Let {L)jjyp and (R)jjyp be such as in Lemma 2.1. There 
exists a constant 9 > such that, for any e,S satisfying 



ed < 9, 



we have 



< 



( L )Hyp + ( R )Hyp 
9 [ x{t)+5e ,(\x-X{t)\ 



9 f*W+ 6s /\ X -X(t) \ , f\x-X(t)\\ ITTI , n , ,, 2 , 
-/ ^ J — K J y - J1 )(U(t,x)-S(t,x)) 2 dx. 

£ JX(t)-8e \ £ J V £ J 



IX(t)-Se 

Proof. We use the definition of X(t) to get 

rX(t) £).4- x ■ Y:/! ^ >' " ! ~ X(l] ' 

:(t 



(L) Hyp = [ X(t) (g)-^( X + X ® )-#( X + X{t) ). V (U\C L ).h(t,x)dx 
yy Jx(t)-Se e \ e / \ e / 



where h(t,x) := [f(u(t, X(t)), - f(U(t, x), C L ) 
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In order to make the function h(t, x) strictly negative over the domain of 
the above integral, we use the condition (j x -Uo)+ € L 2 (R). Indeed, we observe 
that, for any x G [X(t) - Se,X(t)}, 

rX(t) f X(t) 

U(t, X{t)) - U(t, x) = / (8 x U)(t, y)dy < / (d x U)+(t, y)dy 

J x J X 

< \\(d x U(t)) + \\ L 2 m ^\X(t)-x\ < ||(A;7 ) + || L2(R) Vfe, 

where we used that || (d x ll(t))+ 1|^2 is not increasing (see Lemma 3.2). 
We can rewrite the function h as 

h(t,x) = f(u(t,X(t)),^±^) -f{u(t,x),^±^) 
+ f (U(t, x), ° L+ 2 Cr ) f(U(t, x),C L ). 
Since / is increasing with respect to the first variable, we have 
h(t,x) < f(u(t,x) + \\(±U ) + \\ L 2 m VS~e, Cl + Cr ^) - f(u(t,x), ° L+ 2 Cr ) 

+ f (U(t, x), CL+ 2 Cli ) f(U(t, x),C L ). 
Then, thanks to Lemma 3.1, we get 

h(t,x) < A\\(± Uo ) + \\ L , m VS~e- ° L 2A Cr <-0<0 

for \/& small enough. 

Since </>(•), (/)'(■) and r](-\-) > 0, we get 

yy JX{t)-8e £ V £ ' V £ ' 

Then, from Lemma 3.1, we have (changing the constant 8 if necessary) 

yF Jx(t)-8s £ V £ ' V £ ' 



IX(t)-Se 

In a similar way, we obtain the following estimate on (II)h vp - 

r X(t)+Se 
IX(t) 

Combining the two last inequalities gives the desired result. □ 



/ D x ^ f X{t)+5£ ,2.jx-X(t)\^ l( x-X[t)\ . _ . 2 , 

{R) Hyp <-6 J {-)<}>{ — j ±1 )4> f { — j^)(U -C R ) 2 dx. 
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4 Control of the parabolic terms 

For any S > 1/9, we now fix the function <j>$ in the following explicit way. 



(j> s {x) = 



9e 1 ~ ed x, for x e [0, 1/0), 

3 9(x-«) j for xe [i/e,6\. 



(19) 



It can be proved that this functional minimize 



[ (<t>'{x)) 2 XW>e<i>}dx 
Jo 



among all increasing Lipschitz functions with (f>(0) — 0, <f>(6) = 1 and 5 > 1/8. 
We will only use the straightforward computation: 



/ (4>'s(x)) 2 X{4>' s >9<t> s } dx = Cg- e 
Jo 



285 



(20) 



This section is dedicated to the proof of the following proposition. 



Proposition 4.1. There exists constants 9,C > such that for any e,S veri- 
fying 

— < 5 and e5 < 9, 




we have 



dH(t) 
dt 



< Ce 



-68 



Proof. First, we estimate the term (L)j^-j. Integrating by parts, we obtain 

rX[t) '-x + x{t)\,,/-x + x{ty 



2e 



y w- 

J — oo L x 



-X + X(t) 



)d x U( V '(U)-r}'(C L ))dx 
)]\"(U)\d x U\ 2 dx 



Then, using Holder's inequality and Lemma 3.1, we get 



( L )Dif- A 



< 



2e f x ^ 



±±S>)] 2 | 4 I/| ! «ix 



< 



A f x{t) 
8e 

2s r V(M 

A 

C f x{t) 



/^(z£±M) (l/([rW(Ct)) 

y k 

^ — oo 



dx 



-x + X(t) 



\d x U\ 2 dx 



\U -C L \ 2 dx. 
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In the same way, we obtain the following estimate for (Kjjj^j. 



{R) Dlf <- r (— T^J i 

■' e Jx(t) 1 v £ /J 



:(t) 

Combining the two last inequalities, we find 

Q rX(t)+6s 



\U-C R \ 2 dx. 



(L) mf + (R) Dif 



< 



X(t)-Se 



4>'( \- *®\ yf\U{t,x)-S{t,x)\ 2 dx. (21) 



Using Lemma 2.1, Proposition 3.3, and (21), we find 



dH(t) 1 



dt 



L 



X(t)+Se 



e JX(t)-8e 



W - e<t>) 



\x-X{t)\ 



\U{t,x)- S(t,x)\ 2 dx. 



(22) 



Using that U — S is a bounded function, and doing the change of variables 
z = (x — X(t))/s, we find: 



dH(t) ^ C [ XW+5e r 2 

e JX(t)-5e 1 

c\\ u{t)-s{t)\\U f x(t)+S£ 

£ JX{t)-Se 



\x-X(t)\ 



< 



{<t> ) X{C4>'-e<P>o} 



\U(t,x) - S(t,x)\ 2 dx 
\x-X{t)[ 



dx 



<C f ((t>') 2 (z)x{c4>'-e4»a}(z) dz. 
Jo 

Changing the constant 9 if needed, and using (20), gives the desired result. □ 

5 Proof of Theorem 1.2 

This section is devoted to the proof of Theorem 1.2. Integrating the estimate 
of Proposition 4.1 between and t € (0,T) gives the result of (I). Indeed, for 
any e, 5 with | < 5 and s5 < 8, where 9 is the constant from Proposition 4.1, 
we have 

/ r](U(t,x)\S{t,x))dx < H(t) < H(0)+ f ^H(s)ds 

J{\x-X(t)\>5e} Jo dt 

[ v(Uo\So) 



< 



l dx + CTe 



-65 



By taking e := 9 2 , we have for any e < a < e , 



f r](U(t,x)\S(t,x))dx < [ r/(C/ |5o)dx + CTe- Q / £ . 
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It proves our main estimate (9) by taking C* large enough. 
Observe that 



Tj(U\S) dx 



J 

J{\x-X 



(t)\>C*a} 



r]{U\S)dx 



f 

J\\x 



{\x-X{t)\<C'a} 



rj{U\S)dx 



and the second term is bounded by CC*a. Thus, by taking a = elog(l/e), we 
obtain for any t e (0,T), 



r](U\S)dx < [ f](U \S )dx + C*elog(l/e) 
i Jr 

for any e < Eq (changing Eq and C* if needed). 



(23) 



It only remains to prove (11). We define first ip(x) :- 



if > 2, 

1 if |ie| < 1 

2 - \x\ if 1 < |a;| < 2 

Let s e (0,i) and R > 0. We multiply $ R (s,x) := ip C'^ ) to the equation 
(1) and integrate in x to get 

= --^ Jv R -Udx + J d x {^ R )A(U)dx + J d t (* R )Udx + s J * R -dl x Udx 



= -dSl^- 



X(s) 



R 



) • U(s,x)dx 



(i) 

x-X(s) 



R 



) • (a(U(s, x)) - X(s)U(s, x)) dx 



— £ 



R 



(ii) 

^'( X ~* {s) )-d x U(s,x)dx. 



(in) 

By using the above observation, we have 



- X(s)) = Cl ]_ Cr (A{C L ) A(C R ) {C L - C R )X{s)) 



Then we integrate the above equation in time on [0,t] to get: 



\at-X 



(t)\<c( 



+ 



t ■ max 

e(o,t) 



f\l)d t 
Jo 



A(C L ) - A(C R ) - (C L - C R )X(s) - (II) 



(ii') 



(24) 



t ■ max 

se(o,t) 



(^)|)- 



14 



We observe 



(//') < A(C L ) - A(C R ) - 



R 



^ { ^x<A yA{u)dx 



R 



+ 



(C L -C R )X{s) + - I ^( 



„x-X{a). 



R 



)■ (X(s)U{s,x))dx 



For (//(), we compute 



i f -R+X(s) 1 , 

= A(C L ) - / A(U)dx A(C R ) + - / 

ft J-2R+X(s) K Jh 



2R+X(s) 



R+X(s) 



A(U)dx 



< 



1 

R 



r 2R+X(s) 



\A{C L )- A(U)\dx+ / \A{U) - A(C R )\dx 

-2R+X(s) JR+X(s) 



\U-S\dx. 



-R+X(s) 
-2R+X(s) 

We use \A(y) - A(z)\ < C\y - z\ for \y\, \z\ < M Y to get 

(J r2R+X(s 
R J-2R+X(s) 

We use the Holder's inequality and Lemma 3.1 to get 

(H[f<^ - I v(U(s)\S(s))dx. 
ft Jr 

Likewise, for the second term (II^), we have 



1 



(c L -c R ) + - I v'( 



{II> 2 ) = \X(s 

C -\ 

R J-2R+X(s) VR 

where we used |X(s)| < C. Thus we get 



, x-X(s) 



R 



) • U(s,x)dx 



< 



2R+X(s) p 

\U-S\dx< — -\\U( S )-S( S )\\ LHm 

2R+X(s) 



(Il'f<^ - [ v(U(s)\S(s))d2 
ft Jr 



(25) 



On the other hand, we compute 

X(t) 



) • U(t,x)d 



"-hi 



) • Uo(x)dx 



J\l)ds\ = \Ji>(- 

I i>( X ~* {t) ) ■ (U(t, x) - S(t, x)) dx + J i> C~* {t) ) ■ S(t, x)dx 
J • S (x)dx - J • (U (x) - S (x))dx 
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Note that / ipi^^ 1 ) ■ S(t,x)dx = J V(f ) • S {x)dx. Thus, we have 



< 



/ R {t) ) ' ( U ^ X) S ^ x ^) dx \ + | / ^(|) ' - 5o ^)) 



We use the Holder's inequality and Lemma 3.1 to get 



(I)ds 2 <CR([ r](U(t)\S(t))dx+ [ 7](U \S )dx). 

V/TO ./TO ' 



\dx 



(26) 



Also, we have 

(777) 

e 

= R 



fn 



e 
R 

-R+X(s) 



,.x-X(s) 



R 



) ■ d x ll(s,x)dx 



d x U(s,x)dx — \ 

-2R+X(s) Jh 



2R+X(s) 



R+X(s) 



d x ll(s, x)dx 



(27) 



<£. 4 .||[/( s) ||^<^£ 



Finally, by using (23), we combine (25), (26) and (27) with (24) to get, for 
anyt e (0,T), 



Wt-X(t)\ 2 <c(^ + R) ■ (J v (Uo\So)dx + elog(±j) + 



C-e 2 -t 2 
R 2 



Since the above estimate holds for any < R < oo, the estimate (11) follows 
once we take R := t 2 ^ (changing C* if needed). 

Acknowledgment:. The second author was partially funded by the NSF. 
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